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Abstract. We study a Hamiltonian describing a pendulum coupled with several anisochro- 
nous oscillators, devising an asymptotic expansion for the splitting (matrix) associated with a 
homoclinic point. This expansion consists of contributions that are manifestly exponentially 
small in the limit of vanishing hyperbolicity, by a shift-of-contour argument. Hence, we infer 
a similar upper bound on the splitting itself. 



1. Main Concepts and Results 

1.1. Background and history. The study of "separatrix sphtting" dates back to Poincare's 
classic Les Methodes Nouvelles de la Mecanique Celeste [Poi93]. 

Starting with Kolmogorov's 1954 note [Kol54], it was proved in a series of papers over 
a period of twenty years that quasiperiodic motions (invariant tori) are typical for nearly 
integrable Hamiltonians [Mos62, Arn63, Mos66a, Mos66b, Mos67], and that motions which be- 
come quasiperiodic asymptotically in time (stable/unstable manifolds) are stable under small 
perturbations [Mos67, Gra74]. 

Arnold [Arn64] described a mechanism how a chain of such "whiskered" tori could provide 
a way of escape for special trajectories, resulting in instability in the system. (A trajectory 
would typically lie on a torus and therefore stay eternally within a bounded region in phase 
space.) The latter is often called Arnold mechanism and the general idea of instability goes 
by the name Arnold diffusion. It is conjectured in [AA68] that Arnold diffusion due to Arnold 
mechanism is present quite generically, for instance in the three body problem. 

Arnold mechanism is based on Poincare's concept of biasymptotic solutions, discussed in 
the last chapter of [Poi93], that are formed at intersections of whiskers of tori. Following such 
intersections, a trajectory can "diffuse" in a finite time from a neighbourhood of one torus to 
a neighbourhood of another, and so on. 

Chirikov's report [Chi79] is a very nice physical account on Arnold diffusion, while Lochak's 
compendium [Loc99] discusses more recent developments in a readable fashion and is a good 
point to start learning about diffusion. Gelfreich's introduction [GelOl] to splitting of separa- 
trices is excellent, and we recommend it to anyone intending to study the topic. From there 
one should advance to [GLOl], which covers more topics with more details. The extensive 
memoir [LMS03] by Lochak, Marco, and Sauzin is written from the geometric point of view. 
It has a historical flavor, making it interesting and accessible to virtually anyone. 

For the separatrix splitting of the periodically forced pendulum, see [IIMS88,SMH91,DS92, 
EKS93]. In the case of the standard map, an asymptotic expression of the splitting has been 
obtained in [Gel99], as the culmination of a series of works starting with Lazutkin's [Laz03]. 
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For the quasiperiodically forced pendulum, several studies exist. We mention [DGJS97a, 
DGJS97b,DJSG99,DG03]. Especially [Tre97,LMS03,RW00,Sau01,RT04] emphasize geomet- 
rical aspects of separatrix sphtting in the quasiperiodic setting. A lower bound on the splitting 
has been obtained in [GGM99b], where the pendulum is coupled to two rotators — one speed- 
ing up and one slowing down. A generalization to the case of several slow rotators can be 
found in [Pro03]. These results extend [CG94, Gal94, Gen95a]. See also [GGM99a]. 

1.2. The model. We consider the Hamiltonian 

V, /, A) = + cos + _ ^^(^^ ^) (11) 

of a pendulum coupled to d rotators, with € := R/27rZ and / € M the coordinate and 
momentum of the pendulum, and i/; £ T'^ := (S^)"^ and A £ R.'^ the angles and actions of the 
rotators, respectively. The perturbation / is assumed to be real-valued and real-analytic in 
its arguments, and A is a (small) real number, whereas g > 0. This Hamiltonian is sometimes 
called the generalized Arnold model or the Thirring model. It is the prototype of a nearly 
integrable Hamiltonian system close to a simple resonance, as is explained in the introduction 
of [Gen95b]. A review of applications can be found in [Chi79]. 
The equations of motion are 

= /, ^ = A, i = g^sm^ + Xd^f, A = Xd^f. (1.2) 

For the parameter value A = 0, which is addressed as the unperturbed case, the pendulum 
and the rotators decouple. The former then has the separatrix flow : M — > given by 

(j){t) = ^>°(eS*) with ^>°(z) = 4arctanz. 

By elementary trigonometry, the odd function possesses the symmetry property 

^\z) =2tt -<^°{z'^). (1.3) 

In the phase space of the pendulum, the separatrix — given by — separates closed trajectories 
(libration) from open ones (rotation). 

On the other hand, : 1^ — > T'^ is quasiperiodic: the vector uj := A{0) = A[t) in 

iP{t) = 7p{0) + ut (mod 27r) 

is assumed to satisfy for some positive numbers a and the Diophantine condition 

\uj-q\> alql"" for q£Z'^\{0}. (1.4) 

Thus, at the instability point of the pendulum, the flow possesses the invariant tori 

To:={(0,V',/,^) = (0,^,0,a;) \ 9 e T^} 

indexed by lo, with stable and unstable manifolds (Wq and Wg , respectively) coinciding: 

W'o''' = {{^,i^,I,A) = {^\z),e,gzd,^^{z),uj) I [-(»,oo], 0eT'^}. (1.5) 

Remark 1.1. The constant g is the Lyapunov exponent for the unstable fixed point of the 
pendulum motion; in the limit s — oo two nearby initial angles (/>(s) and (j){s + 6s) separate 
at the exponential rate e^'*. As (p{t) = <I>'^(e*/^ ), this fixes a natural time scale of g~^ units, 
characteristic of the pendulum motion in the unperturbed Hamiltonian system (II. ip . 

When the perturbation is switched on (A 7^ 0), it is known that most of the invariant tori 
survive and have stable and unstable manifolds — or "whiskers" as Arnold has called them — 
that may not coincide anymore. We prove a bound on their splitting. 
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1.3. Main theorem. In [Ste07] we construct the perturbed manifolds in a form similar to 
(jl.Sp as graphs of analytic functions over a piece of [—00, 00] x T'^. To that end, we look for 
solutions of the form 

(0(t),V'(t)) = ($(e^*,cjt),wt + ^(e^*,cjt)) = (0,wt) + ($,^')(e^*,u;t) (1.6) 

with quasiperiodic behavior in one of the two limits t — > ±00. Note especially that the 
Lyapunov exponent 7 > depends on A, with 7|;^=q = g. 

Remark 1.2. One should not expect asymptotic quasiperiodicity in both limits t — > ±00, as 
the unstable and stable manifolds, and VV|, are generically expected to depart for nonzero 
values of the perturbation parameter A. Therefore, either the past or future asymptotic of a 
trajectory will evolve so as to ultimately reach the deformed invariant torus T\. 

Let us denote the total derivative d/dt by dt and the complete angular gradient (5(^,3^) by 
d for short. Substituting ()1.6p into the equations of motion, we get the equation 

{uj-de+ je^'d,fX{e^\Lot) = [{g^ sin 0) + A df{X + (0, ^))](e^*, 

for X := (<I>, ^'), where 9 stands for the canonical projection [—00, 00] x 
Notice that the partial differential operator 

C := oj ■ 80 + jzdz 

satisfies the characteristic identity 

CF{ze^\ e + iot) = dtF{ze^\ O + ut), (1.7) 
which reflects the time derivative nature of £. In fact, if T is the "time-reversal map" 

T{z,e)^{z-\-9), (1.8) 

then, by the chain rule, 

C{FoT) = -{CF)oT. (1.9) 

Let us abbreviate 

Vt{X):=[g'^sm^,Q) + \9.{X) with 0(X) := + (0, 6*)). (1.10) 

We have then encoded the equations of motion into the PDE 

C^X = VL{X). (1.11) 

The action variables, or momenta, trivially follow from the knowledge of X{z^6): 

{I{t),A{t)) = (0, Lo) + Y{e^\iot), Y :=CX. 

The solutions X provide a parametrization of the deformed tori and their stable and unstable 
manifolds. As hinted below p.6p . we find two kinds of solutions, X'^{z,6) defined for z € 
[— zc^o] and X^{z,9) defined for z S [—oo,—Zq^] U [2;(5~^,oo]. Here, zq > 1. 
We will consider A small in a g-dependent fashion, taking 

e := A5-2 (1.12) 

small. Such a choice is needed for studying the limit g ^ 00, which corresponds to rapid 
forcing; see Remark 11.11 The domain we restrict ourselves to is 

D := {{e,g) E C x M | |e| < eo, < 5 < go}, (1.13) 

for some positive values of eo and go. 

The following theorem is from [Ste07]. It is a version of a classical result, and by no means 
new; earlier treatments include [Mel63, Mos67, Gra74,Eli94, Gal94, Gen95a, Gen95b]. 
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Theorem 1 (Tori and their whiskers). Let f be real-analytic and even, i.e., 

/(</), V') = /(-</), -V). 

Also, suppose oj satisfies the Diophantine condition (jl.4|) . and fix go > 0. Then there exist a 
positive number eo and a function 'y{e,g) on D, analytic in e with \^ — g\ < Cg\e\, such that 
equation (jl.lip has a solution which is analytic in e as well as in {z, 9) in a neighbourhood 
of [-1, 1] X and which satisfies (X^ := ($°,0); 

X-(1,0) = (^,0), X-{z,9)=X\z) + 0{e). (1.14) 

Corresponding to the same 7, there exists a solution X'^{z,9) = X^{z) + 0{e) which is an 
analytic function of {z"^, —9) in a neighbourhood of [—1, 1] x T"^. The maps 

T4/^.«(z,^) = (X".«,y^.«)(z,6/) + ((0,6'),(0,a;)), y^'" := (1.15) 

provide analytic parametrizations of the stable and unstable manifolds W^'" of the torus T\ . 

Remark 1.3. The number eo above depends on the Diophantine exponent i' and on /. The 
perturbation {4>,il^) 1-^ f{(f>,ip) is analytic on the compact set S"*^ x T'^. By Abel's Lemma 
(multivariate power series converge on polydisks), it extends to an analytic map on a "strip" 
|9m(/)|, |3'm'i/'| < (^ > 0) around x T'*. By Theorem[Tl there exists some < a < rj such 
that each 9 1— > X^'^{ ■ ,9) is analytic on |9m6'| < a. 

An important part of Theorem [1] is that the domains of X" and X'^ overlap. Namely, if X 
solves equation (jl.lip . then so does (27r, 0) — X oT. This is due to (jl.9p and the parity of /. 
Consequently, by time-reversal, the stable and unstable manifolds are related through 

X' = {27r,0) - X"" oT. (1.16) 

In particular, as r(l,0) = (1,0), ^''(1,0) = X"(1,0). The actions y*'" = CX''"" satisfy 
ys — y« Q j^^ yielding y'^(l,0) = y"(l,0). In other words, a homoclinic intersection of the 
stable and the unstable manifolds W^'" occurs at {z,9) = (1,0), as their parametrizations 
p.lSp coincide at this homoclinic point. Since the manifolds W^'" are invariant, there in fact 
exists a homoclinic trajectory on which the parametrizations agree: 

W'{e^\ujt) = Wie^^uot). (1.17) 

Coming to the second result of [Ste07], let us expand X" = YlT^o e^X"'^. It turns out that 
the common analyticity domain of each X^'^ in the z-variable is in fact much larger than the 
(small) neighbourhood of [—1, 1] — the corresponding domain of X^ according to Theorem [TJ 
Namely, it includes the wedgelike region 

U^,,?:={kl <t] |J{argzG [-^9, i?] U [vr - i9, vr + CC (1.18) 

(with some positive r and "!?). We repeat parts of the argument in Section [5] for convenience. 

Theorem 2 (Analytic continuation). Each order X"^'^ of the solution extends analytically to a 
common region U,-^^ x {|Qm^| < a}. Moreover, if ip ^ f{ - ^ij:) is a trigonometric polynomial 
of degree N, i.e., N is the smallest integer such that f{-,q) = whenever \q\ > N, then 
9 I— > X"'^( ■ ,9) is a trigonometric polynomial of degree IN , at most. 

Remark 1.4. With r] and a as in Remark 11.31 the numbers r and ?? are specified by the 
following observation: (^^{z) = 4arctan^ implies that |9m$'^(2:)| < rj \n V^-.-d with r and 
sufficiently small. By Remark 11.31 (-2,^) ^ f{^^{z),9) is analytic on U,-^^ x {|Qm0| < a}, 
which is the basis of the proof. 
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In spite of Theorem [21 (a straightforward upper bound on) X"'^ grows without a hmit as 
\^ez\ oo, such that there is no reason whatsoever to expect absolute convergence of the 
series Y^£^o e^X"'^ in an unbounded z-domain with a fixed e. In fact, it is known that the 
behavior of the unstable manifold gets extremely complicated for large values of z even with 
innocent looking Hamiltonian systems. Still, it seems to us that the possibility of a uniform 
analytic extension of the coefficients X"'^ has not been appreciated in the literature. 

Due to (jl.l6p . an analog of Theorem [2] is seen to hold for , with z replaced by z~^. 

Theorem [2] allows one (at each order in e) to track trajectories t i-^ W^''^{e^^ ,9 + tot) on 
the invariant manifolds W^'" for arbitrarily long times in a uniform complex neighbourhood 
|9mt| < g~^^ of the real line, for arbitrary 6 € T"^. The motivation for doing this stems from 
studying the splitting of the manifolds W^'^ in the vicinity of the homoclinic trajectory (|1.17p . 
The general ideology that, being able to extend "splitting related functions" to a large complex 
domain yields good estimates, is due to Lazutkin [Laz03], as is emphasized in [LMS03]. 

In order to study the intersection more closely, we express the actions as functions of 
the original angle variables (</>, ^) = X^'^{z,9) + {0,0) appearing in the Hamiltonian (jl.ip . 
To this end, let F**'" : (z, 6) (</>, ip) be the above coordinate transformations, and write 
Ys,u ^ y^,« o ps,u^ set 5 := {d^, d^) and D := {d-^,de). By the chain rule, 

^ys,n ^ ^ ps,u^^ DF''"". (1.19) 

The invertibility of the matrix DF^''^ is a consequence of Theorem [TJ At the homoclinic point 
{z,e) = (1,0), equation (fTTTH implies that F' {el\ ujt) = F"(e'^*,u;t) and F^'"(1,0) = (7r,0). 
Casting in the obvious manner y"'* = (Y^'^ ,Y^'^), we define the splitting vector: 

A{^, i;) := {A^ - AY = in - nVi^, i^)- 

Transposition indicates that we consider A a column vector. We define the splitting matrices: 

T{t):=de{Y^-Y^f{e'^\u:t), (1.20) 

T{t) := d^{Y^ - nf {F^^'{e'<\Lot)). (1.21) 

Notice that if {(j),ij)){t) stands for the moment for the angles on the homoclinic trajectory at 
time t, with the particular initial condition {(j),tp){{)) = (7r,0), then T = d^A{(j),tp). 

It is nontrivial but straightforward to establish that T and T differ by a close-to-identity 
factor, if e is taken to be small proportionally to g: 

Theorem 3. Fix t G M and define e := g~^e^^^^e. Then, as g ^ 0, 

T(t) = T(t) (l + 0(e)), 

if e is sufficiently small (independently of g and t). 

Theorem [3] is our first result. Its proof in Appendix lAl is based on energy conservation and 
the fact that the actions are given by scalar potentials. We can now state the main theorem: 

Theorem 4 (Homoclinic splitting). Let us assume that tl^ ^ f{ - is a trigonometric poly- 
nomial in addition to the assumptions of Theorem [21 Then for each f € M there exist positive 
constants C and c, such that the exponentially small upper bound 

\T,,{t)\<C\e\e-^'-'''^^\ 
where v is the Diophantine exponent, holds, provided eg~^ is sufficiently small. 
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Theorem m is derived from Proposition 13.121 — the centerpiece of this work — stated in Sub- 
section [3^ Technicahties aside, we beheve that Subsection l3.3l can be understood without any 
further introduction, and we urge the reader to go through it before moving on to Section [2l 

Why just study a d x d sub matrix of the fuh (d+l) x (d+l) Jacobi matrix of ^4*) 
when measuring transversality of the homochnic intersection PI W]^? Because the latter is 
singular. As a matter of fact, ()A.2h simply means that the splitting vanishes in the direction 
of the homoclinic trajectory. The author is grateful to Dr Mischa Rudnev for pointing this 
out. For further motivation, see item Rl of Appendix R and Section 8 in [GGM99b]. 

Acknowledgements. I am indebted to Antti Kupiainen for all his advice. I express my grat- 
itude to Guido Gentile, Kari Astala, and Jean Bricmont for their comments. I thank Giovanni 
Gallavotti and Emiliano De Simone for discussions at Rutgers University and University of 
Helsinki, respectively. Mischa Rudnev, Pierre Lochak, and Michela Procesi were kind enough 
to explain me some of their works on the subject. I am grateful to Joel Lebowitz and Rutgers 
University for hospitality during the final stages of writing this work. This work was supported 
by the Finnish Cultural Foundation and NSF Grant DMR-01-279-26. 

2. Analytic Continuation of the Solution 

Here we present the proof of Theorem [2l We drop the superscript u from the notation and 
single out the uncoupled part of the complete solution X; 

X = X° + X with Xl^^ = 0. 
Equation (fTTT]) and C^X^ = {j'^sm^^,0) imply C'^X = -(7^ sin ^>o, 0) + n{X^ + X). Thus, 

ICX = W{X), (2.1) 

when we define the linear operator 

IC:=(^ ^2) with L := £2 _^2^os$° (2.2) 



.0 ^ . 
and the nonlinear operator 

W{X) := (-7^ sin $° - 72(cos $°)$, 0) + 17(A:° + X). (2.3) 

Throughout the rest of the work, we shall refer to different parts of the Taylor expansion 
of a suitable function h{z, 9) around z = using the notation 



hk{9) := h^k{z,e):=Y,^''f^k{0) and d^h := h - h< 



kl 

By dZl]), 62X satisfies 
where 



j=0 



IC62X = Wi62X), (2.4) 



W{Z) :=62 



(2.5) 



Let us consider analytic functions Z on a small, compact and complex, neighbourhood H 
of the set [-1, 1] x T'^. It is a complex Banach space A, once equipped with the supremum 
norm, and has the closed subspace 

Ai:={Z eA\Z<i=0}. (2.6) 

We have shown in [Ste07] that /C maps ^1 to itself and has a bounded, 0(7"^), inverse. 
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By virtue of ()2.4p and W^s analyticity, 62X admits the representation 



d2X = IC-H2 



7^cos<I>° 0\ ~ rfcW^ ^ 



fc=0 



+ 



00 

+ /C"^ > \w 

k=l 



00 



(2.7) 



on the set H, taking e small enough, and denoting 

1 



w 



A;! 



(2. 



as well as a repeated argument of such a symmetric fc-linear operator by (x)^'' := {x,. . . ,x), 
for the sake of brevity. Observe that we have omitted a 82 in front of the square brackets on 
the second line of (|2.7|) as redundant. 

Equation ()2.7p may be viewed as a recursion relation for 62X. It is crucial that 



(2S 



when {e,g) G D; see ()1.13p . Namely, any given order ^2-^ in the convergent expansion 62X = 
Yl'^i £^ ^2X^ is then completely determined by X<i and the lower orders 62X'' {1 < I < i — 
through the right-hand side of (12. 7p . Moreover, since X<i = 0(e), only finitely many terms 
in the sum over the index k are involved. Together these facts imply that only finitely many 
recursive steps using (j2.7p are needed to completely describe any given order S2X''' in terms of 
X<i alone and that, at each such step, only finitely many terms from the fc-sum contribute. 

It is important to understand that X<i is a predetermined function. As we shall see, the 
recursion procedure will then provide the analytic continuation of each X"'^ = X^^ + 62X^ 
(^ > 1) to the large region Ur,,? x {|9m^| < a} of Theorem [2l 

2.1. Tree expansion. We next give a pictorial representation of the above recursion. It 
involves tree diagrams similar to those of Gallavotti, et al. (see, e.g., [Gal94, CG94]), with 
one difference: there will be no resummations nor cancellations, as the expansion in p.7p 
contains no resonances and is instead well converging. This so-called tree expansion is needed 
for bookkeeping and pedagogical purposes; we simply choose to draw a tree instead of spelling 
out a formula. 

Let us first define the auxiliary functions 



w 



(0) 



, (g>fc 



if /c = 2,3, 



and make the identifications 

:= IC-^62h^'''^ and 



(2.10) 



A:=0 



Furthermore, let 



-O -=^2^' 



X 



<1) 



and 
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In the diagram representing the A:-hnear w^^\ the k "free" Unes to the right of the node stand 
for the arguments. We say that these hnes enter the internal node, whereas the single hne to 
the left of the node leaves it. For instance, 




X<i, 62X, JC~^62h^^^). 



Notice that, as w^^'^ is symmetric, permuting the lines entering a node does not change the 
resulting function. We emphasize that all of the functions introduced above are analytic on 
n and |e| < eo. 

In terms of such tree diagrams, or simply trees, equation ()2.7p reads 




(2.11) 



using multilinearity to split the sums X<i + 62X into pieces. Above, the sum after the first 
tree consists of all trees having one internal node and an arbitrary number of end nodes, at 
least one of which, however, is a white circle. This rule encodes the fact that on the second 
line of ()2.7p the summation starts from k = 1 and that the contributions with only X<i in 
the argument {i.e., trees with only black dots as end nodes) are cancelled. 

Using ()2.7p recursively now amounts to replacing each of the lines with a white-circled end 
node by the complete expansion of such a tree above. This is to be understood additively, so 
that replacing one end node, together with the line leaving it, by a sum of two trees results in 
a sum of two new trees. For example, such a replacement in the third tree on the right-hand 
side of (|2.1ip by the first two trees gives the sum 

Before proceeding, we introduce a little bit of terminology. The leftmost line in a tree is 
called the root line, whereas the node it leaves {i.e., the uniquely defined leftmost node) is 
called the root. A line leaving a node v and entering a node v' can always be interpreted as the 
root line of a subtree, the maximal tree consisting of lines and nodes in the original tree with 
V as its root. We call v a (not necessarily unique) successor of v' , whereas v' is the unique 
predecessor of v. 

The recursion ()2.1ip can be repeated on a given tree if it has at least one white circle left. 
Otherwise, the tree in question must satisfy 

(Rl') The tree has only filled circles (^^) and black dots (•) as its end nodes, 
together with 

(R2') Any internal node has an entering (line that is the root line of a) subtree containing 
at least one filled circle as an end node. 
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After all, the recursion can only stop by replacing an existing white circle with a filled one. 
Continuing ad infinitum yields the expansion 

Q = (Trees satisfying (Rl') and (R2')) = T, (2.12) 



trees T 

where the prime restricts the summation to trees T satisfying (Rl') and (R2'). We point out 
that each admissible tree appears precisely once in this sum, considering different two trees 
that can be superposed by a (nontrivial) permutation of subtrees that enter the same node. 

The earlier discussion concerning the description of 52X^ in terms of a finite sum involving 
only X<i translates to the language of trees in a straightforward fashion. First, the second 
part of (|2.9p and X<i = 0{e) amount pictorially to 



0(e) and 



Oie), 



since /C ^ produces a factor of g ^. Second, w^''^ = 0{g'^) and the first part of (j2.9p yield 



= Oie') {k > 1) 



and 




0(1) (k > 2). 



Expanding the filled end nodes 



oo 

E 

k=l 



(2.13) 



according to (j2.1Up . on the right-hand side of (|2.12p . we get a new version of the latter by 
replacing the rules (Rl') and (R2'), respectively, with 

(Rl) The tree has only numbered circles ((^ with arbitrary values of k) and black dots (•) 
as its end nodes, 

and 

(R2) Any internal node has an entering (line that is the root line of a) subtree containing 
at least one numbered circle as an end node. 

Let us define the degree of a tree as the positive integer 

degr:=#( . ) + #( .) + 



oo 

E 

k=l 



(2.14) 



for any tree T satisfying (Rl) and (R2). By t^(G) we mean the number of occurrences of 
the graph G in the tree T. That is, the degree of a tree is the number of its end nodes with 
suitable weights plus the number of nodes with precisely one entering line. Since a tree has 
finitely many nodes, its degree is well-defined. Then a rearrangement of the sum arising from 
(|2.13p being inserted into (j2.12p yields formally 



oo 

E 

1=1 trees T 
deg T=l 



E' 



(2.15) 



where the asterisk reminds us that the rules (Rl) and (R2) are being respected. 

According to the analysis above, the particular graphs appearing in the definition of deg T 
are the only possible single-node subgraphs of T proportional to a positive power of e. Since 
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each tree is analytic in e, writing for the A;th coefficient of the power series in e, we have 

oo oo 
A;=deg T fc=0 

Hence, only trees with degree at most equal to i can contribute to 62X^: 



= E E = ( E ^) (2-16) 

l=ldegT=l degT<i 



or, alternatively, 



52X= r + 0(e^+i) (e^O) (2.17) 

deg T<1 

for each and every i = 1,2,.... The expansion in (|2.15p is just a compact way of writing 
(|2.17p . We emphasize that the latter can be derived completely rigorously, for each value of 
i separately, but resorting to the use of formal series allowed us to treat all orders of at 
once. We call the series (j2.15p an asymptotic expansion of 62X; the partial sums YldcgT<e'^ 
need not converge to 62X for any fixed e as i ^ 00, but for a fixed i the error is bounded by 
an ^-dependent constant times |e|^^^ on the mutual domain of analyticity, |e| < eo- 

Example 2.1. The beginning of the asymptotic expansion ()2.17p reads 

62X = Q+0{e') = + + . + 




+ 0(e3). 

2.2. Analyticity domain of trees. As already pointed out, all trees T above are analytic 
functions of {z,9,e) on 11 x {|e| < eo}. Due to the projections 62 appearing in ()2.10p . they 
also satisfy T\z=q = dzT\z=o = 0, i.e., are elements of the space ^1 defined in ()2.6p . On this 
space, the operator IC in ()2.2p has a bounded inverse given by a simple integral kernel [Ste07]. 
Consequently, the analyticity domain of a tree in the z-variable is in fact much larger than 
the neighbourhood of [—1, 1] included in 11; it contains the wedgelike region Vr.^ in (ll.lSp . 

Lemma 2.2 (Analytic continuation of trees). Without affecting the analyticity domain with 
respect to e, there exist numbers 0<r< 1, 0<'i?<7r/2, and < a < r] such that each tree in 
the sums (j2.15p and (|2.17p extends to an analytic function of {z,9) on U,-^^ x {IQ'm^l < o"}. 

Since the sum in ()2.16p is finite and the functions X<i and X^ in A = A'^ + A<i + 62X 
are analytic on Ut-,i9 x {|9m0j < a}, Theorem [2] follows from Lemma 12.21 

3. Size of the Homoclinic Splitting 

The Lyapunov exponent 7 = j{€,g) is an analytic function of e in a neighbourhood of the 
origin (|e| < eo). In fact, 7 = g{l + 0(e)) for small e. Further, since C depends on 7, we have 
C = Co + 0{e) with = u! ■ de + gzdz ■ In particular, as X^{z,6) = (4 arctan z, 0)'^, the 
expansion of the splitting matrix T — see ()1.20p — begins linearly with respect to e: 

00 

T = ^e^T^. 

1=1 

It is customary to call the first order coefficient, T^, the Melnikov term or the Melnikov matrix. 
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3.1. The Melnikov term. Let us study T^. The symmetry (|1.3p comes in very handy, since 
it imphes, together with the even parity of /, that f{^^{z~^), —9) = f{^^{z),6). Setting 

/oo 
a2/($0(e5*+«), e + ivt) - dlfiO, e + ivt) dt (3.1) 
-oo 

then yields, after some work, 

T\t)=g''F{0,0) 

for all t. Here the second term in the integrand removes the quasiperiodic limit of the first 
one, making the integral absolutely convergent. In fact, the integrand tends to zero at an 
exponential rate as |t| — oo. 

We now give a result, due to Lazutkin in it's original form, [Laz03], and extended by 
[DGJS97a], [SauOl], and [LMS03] to the quasiperiodic setting. See Appendix 1X1 for the proof. 

Lemma 3.1 (Lazutkin). Suppose a function F is analytic on S := [—i-d,i'd] x {|9m0I < r]}, 
continuous on the closure S, and satisfies the identity 

F{s,9) = F{0,6 - g-^LUs). (3.2) 

Then F extends analytically to {\Qms\ < -i?} x {|9m0| < rj}, with (|3.2|) still holding. On 
M X T*^ obeys a bound of the form 

\F{s,e)-F\<CB{g)e-'^^-''"'^'\ B{g) := sup \F{s,e)l 

is,e)es 

where F stands for the 9-average {F{s, ■)) and is independent of s. 

Remark 3.2. Notice that, for small values of g, the analyticity domain of F in Lemma l3. II is 
much larger than what the right-hand side of (13. 2p suggests. 

The lemma applies to the function defined in (|3.ip . due to Remarks 11.31 and 11.41 regarding 
the analyticity of /. Because the ^-derivatives in the integrand are interchangeable with total 
^-derivatives, F = 0, and we get 

Proposition 3.3. There exist positive constants c and C such that the Melnikov term satisfies 

\T'\<Cg^e-'^<^-''"'^'\ 
Remark 3.4. Identity (13. 2p in Lazutkin's lemma is equivalent to 

ids + g-^u ■ de)F = 0, (3.3) 

or, calling F{z, 6) = F{g~^ In z, 9), to CqF = 0. There exists a whole industry trying to push 
through the argument given above replacing with a "better measure" of the splitting, such 
as T, by searching for a coordinate system in which the "measure" satisfies ()3.3p . The state 
of this quest is best described in [LMS03]. 

Neither Lemma l3.1l nor Proposition 13.31 is optimal. There is a refinement in [SauOl] which, 
however, is not optimal. In many special cases one can derive a much stronger bound. 

3.2. Regularized integrals. This subsection borrows heavily from Gallavotti, et al; see for 
instance [CG94] and [GGM99b]. We adhere to the notation 

X{t;z,9) := X {ze^\ 9 + ujt) , 

which makes the dependence on the Lyapunov exponent 7 implicit (Recall that 7 is a function 
of e), such that the perturbation expansion reads compactly X{t;z,9) = Yl'^Qe^X^{t; z,9). 
With little trickery, the equations of motion retain their appearance: defining i}(X)(t; z, 9) = 
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[{g^ sm^,0) + Xn{X)]{t; z,9) together with n{X){t;z,e) = df{X{t;z,e) + {0,9 + cut)), and 
implementing ()1.7p . we obtain the analogue of (ll.llj) : 

dtY{t;z,e) ■.= d'jX{t;z,e) = n{X){t; z,e). (3.4) 

There is a need to consider functions h{t; z, 9) that can be expanded as finite sums 

p 

^tPhP{ze^\d + ujt), (3.5) 

p=0 

where is analytic on {0 < \z\ < r} x T'^ — a punctured neighbourhood of the origin times 
the d-torus — with at worst a finite pole at z = 0. 

Let us set TNh{t; z,6) := J2\q\<N ^(^'i ^)^^'^^ ^^'^ define a regularized integral 



/ 



1 /■* 

h(t;z,e):= lim res-/ e~-^l^lr7v/i(r; z, 6*) dr, 

Af^oo _R=0 R 



the residue at i? = meaning that of the analytic extension from the complex half-plane 
with 3f?c R sufficiently large. This is well-known from the theory of Laplace transforms. The 
truncation, T^, is needed to insure that the origin of the -R-plane is not an accumulation point 
for poles, since factors {R + iuj ■ q)~^ appear from the integral. 

By the following proposition, it is natural to employ the compelling notation 

h{T;z,9) dr := ih{t;z,9), 

— oo J 

and to view z and 9 as mere parameters: 

Proposition 3.5. The regularized integral is linear. Moreover, 

(1) f h{t; z,9) = f h (to; z, 6) + J^* /i(t; z, 6) dr for any real to, 

(2) if h and h' are trigonometric polynomials (in 9) and h{t;z,9) = h'{t;z,9) for all t, 
then f h {t; z,9) = f h' {t; z, 9) for all t, 

{3)Tp^fh = f TNh, 

(A) dtfh{t;z,9) = fj'_^h{T;z,9)dT = h(t;z,9), 

(5) ff h {t; z,9)=: £^ h{s; z, 9) ds dr = £^(t - r) h{T; z, 9) dr. 

Proof. For a finite to, resij=o-^ /j* e~^l'^lti(T) dr equals the evaluation J^* e^^l'^lii(T) c^t|^^q = 

J^* u{t) dr, since the integral is entire in R, from which (1) follows. If is a trigonometric 

polynomial, then f h {t; z, 9) = res j-^a ^^^'^'/^(t; z, 9) dr, yielding (2). Claim (3) is trivial 

and (4) follows from (1). In order to prove (5), one first checks that it holds for trigonometric 
polynomials and then deals with the general case. □ 



We also define 



/ h{T;z,9)dT := lim res-/ e'^^^^TNhiT; z,9) dr 

Jt N^oo R=0 R Jf 



when it makes sense, and write 



ft /'itoO ft fOO POD 

f =-f I +1 =1 . 

J±oo Jt J —oo Jt J —oo 
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The operator D is defined by giving the expression of its Fourier kernel: 

'(iio-q)-^ ifp = gEZ'^\{0}, 
otherwise. 



V-\p,q) 

It is the formal inverse of T> on functions h{z, 6) with vanishing ^-average. Finally, 

.0 

Ih{z,9):= / /i(ze^^ ,9 + lot) dr. 



Applying ()1.7p in the case F = Ih, the formal identity CI = 1 follows. Recalling 6ih{z,6) :- 
h{z, 6) — ho{6) = h{z, 9) — h{0, 9), one also gets IC5i = (5i by a direct computation. 



Proposition 3.6. Let h be analytic on {\z\ < r} x T^, and denote h{t;z,9) = h{ze^* , 9 + uot) . 
ft 

h{T- z, 9) dr = {ho)t + V-^ho{9 + ujt) + 16ih[ze^\ 9 + ujt) 



f 

J — ( 



-oo 

and 



f 

J — ( 



drh{T; z, 6) dr = h{t; z, 0) — {ho)- 



Proof. Decompose h = (ho) + [/iq — (/iq)] + 6ih and use the linearity of For the second 
identity, notice drh{T; z,0) = 'Dho{0 + ujt) + C5ih{ze^'^ ,6 + lot), and recall the identities 
p-ip/io = ho- {ho) and IC6ih = Sih. □ 

Remark 3.7. Observe that, unless the integrals converge in the traditional sense, f'i^ does 
generically not tend to £^ or as t tends to oo or — oo, respectively. 

A detailed proof of the following elementary result can be found in [Ste06] . 

Proposition 3.8. Let So and Si be the spaces of analytic functions h on {\z\ < r} x T*^ 
satisfying {ho) = and {hi) = 0, respectively. Denote h{t;z,6) = h{ze'^^ , 9 + cot) . The 
operators Lt = — 7^ cos ^^{ze"*^), L = £^ — 7^ cos , and C = T) + jzdz satisfy 

dth{t;z,9) = Ch{ze^\9 + ujt) and Lth{t; z , 9) = Lh{ze^\ 9 + ujt) . 

The operator f maps So into itself and, in fact, ^ = d^^ on So- Similarly, the operator 
defined by K^{t, r; z) h{T; z, 6) dT maps Si into itself and is the inverse of Lt on Si. In 
particular, the operator ^ = ( g ^2 ) = ( ^* 92 ) is invertible on the space Si x So, where 

lC-^h{t; z,e) = -f K{t, r; z)/i(r; z, 9) dT. 

Above, the kernel K = ( ''g* ) splits into the useful sum 

1 2 

K{t,T-z) = Y,Y.^t - TYK,j{t-z)Ki,{T-z), (3.6) 

i=0 i=l 

setting P{t]z) := (z^e^T* + l)~^ze'^\ Q{t;z) := z'^ e'^^z"^ e^^^ - 1), and 

^°i = ^(o 0)' Koi = Kii={^^ Ko2 = -^(q 
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3.3. Asymptotic expansion for the splitting matrix. Starting from (j3.4p and using 
Proposition [3S] with {Y(^) = (VX^) = 0, we liave 

z,e) = -f fi(X")(T; z, B) dr 

The identity = (2^,0) - X'^i-t; z-\ -6) yields Y'{t;z,e) = ¥""{-1; z-\ -6), such 

that defining the d component column vector A(t; z,9) := {Y^ — Y^){t; z, 9) and the functions 
f^''^{t; z, 6) := z, 9) + (0, 9 + ojt)) thus results in an expression for the splitting matrix 

T of ()1.20p in terms of regularized integrals: we get 

T(t) = aeA(t;l,0) = de^{Q\e'\ujt) (3.7) 

because 

/t peso 
dgf;{T;z,9)dT + Xf dgf;iT;z,9)dT, 
-oo J t 

where the subindices attached to stand for partial derivatives (our convention is cp = (</>, ip) 
and d = d^ = {d^,d^)). 

Remark 3.9. Extracting here the first order in e = g~^X casts the Melnikov matrix in a 
compact form in terms of the regularized integrals: 

/oo 
dlf{<!>'{en,u^r)dT. 



oo 



The reader is invited to compare this with (|3.1|) and the equation following it. We point out 
that a similar procedure, using Lemma [3. 11 that resulted in the exponentially small bound on 
can be applied here, despite the unusual integrals. 



Suppose that, even with the superscript u, the integral 



/oo 
defl{T;z,9)dT (3.8) 



Jt 

exists. If the integrand is a trigonometric polynomial in 9, then Proposition 13.51 implies that 
the value of the intcCTal only depends on the integrand's restriction to {{t;z,9) \ t € M}, as 
the notation suggests □. Now, fixing {z,9) = (1,0) and dropping it from the notation, 

/oo f'OO 
def:^{T)dT + Xj- [fl,^de{X'-X^)]{T)dT, (3.9) 
"OO J t 

since X"(r;l,0) =X'^(r;l,0). 

Let us make the assumption that ip t-^ f { ■ , ip) is a trigonometric polynomial, such that, by 
Theorem [21 each order of X in e is a trigonometric polynomial in 9. Even so, equation (|3.9|) 
is formal in the sense that the integrands are defined for r ^ only up to an arbitrarily 
high order in e. However, it is asymptotic, which is to say that at each order the identity is 
exact. Put differently, (|3.9|) is a collection of exact identities, one for each order in e, written 
in closed form. Moreover, at each order, the integrands are analytic functions by virtue of the 
extension result in Theorem [2j Of course, we need to check that each order of ()3.8p . i.e., the 

integral [def^Yir; z, 9) dr with arbitrary ^ G N, exists 0. 

The point of the formula above is twofold. First, the integral £^ dgf^{T) dr (at each order 
in e) turns out to be exponentially small in the limit (7 — > 0. Second, we can actually construct 



^This is the first of tlie two places wiiere we need to restrict ourselves to trigonometric polynomials. . . 
^. . . and this is the second. 
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a (formal) recursion relation for the function dQ{X^ — X")(r), taking us to ever-increasing 
orders in e, as follows. Differentiating both sides of ()3.4p with respect to one obtains 

with Lt = df — 7^ cos <I>''(2;e'^*), as in Proposition 13.81 ^^'^ 

Owing to the ^-derivative acting on X^'" and Proposition 13.81 we get 

deX^'^it; z, e) = f K{t, t'; z) [Af ^'"^e^^'" + \Q] {t'; z, 9) dt' , 

upon choosing oo in conjunction with X'^ and —oo with X^. The most convenient way to 
see this is to infer the identity for dgX^ and employing K{—t,—T;z^^) = —K(t,T;z) and 
X'it; z, 9) = {2tt, 0) - z-\ -9). 

At {z,9) = (1,0), M'{t) = Wit) and f^^{t) = f^^{t), such that (2) of Proposition [33] 
validates 



r — oo 

de{X^ - X")(r) = / K{t, t') [M^deX^ + Xf;^] (r') dr' + 

J oo 

+ f K{t, t') [M^deiX' - X-)] (r') dr' . 

J oo 



(3.10) 



By construction, the quantities appearing in square brackets above are ^-gradients: 

M^'-deX^'- + A/^; = de^''^ 

with, for instanced, 

o,s,u . (Lt 0\ ^ „ c,(Y'^^u\ /^7^cos$°(ze'^*) 0\ ^ „ 

^ -(,0 d^)^ -^^^^ >~\ oj^ ' 

regardless of the value of [z,9). In addition, but only at {z^9) = (1,0), 

[M'de{X' - X")] (t) = [9,(21^ - 21")] (t). 

Notice that ()3.10p is a fixed point equation of the form C, = £ + BM^( solved formally by 
C(t) = dg{X'^ — X^){t). This is highly interesting from the point of view of asymptotic analysis, 
since = 0(e) multiplies C on the right-hand side. Indeed, truncating the Taylor series in 
e, e.g., C-^ = E(Lo the identities ('' = {£ + BM%f = {£ + BM%^^-^f become exact. 
Hence, we can recursively construct any order of dQ{X^ — X'^){t) from its lower orders, without 
the need of diverting from {{t\z,9) = (t;l,0) | i G M} in order to compute the ^-derivative. 
Cumulatively, we have C,^ = [X^^Zq (-BM'^)-'f ] where j is a power. 

Proposition 3.10. In brief, the (presumably) divergent Neumann series '^J^^iBM^y £ is 
an asymptotic expansion of dQ{X'' — X^)[t) in the sense that 

(1) At each order in e it terminates after finitely many well-defined terms, and 

(2) \de{X' - X") - [YT=Q{BMy£] (t) = 0{e''+^) for each A; G N. 



Once inserted into ()3.9p . the latter series provides us with an asymptotic expansion for the 
splitting matrix T in ([31 



"^We could make 21°'" of order e by deleting its lowest order, because the latter does not depend on 9. 
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Corollary 3.11. In the asymptotic sense of Proposition [3.10\. 



T{t) = xj dgf^ir) dr + xi f^^ Y^^BMyS (r) dr. 

J -co Jt '- j_Q 



The operator B above has the expression 

Bh{t;z,e) = S K{t,T;z)h{T;z,6)dT, 

J oo 



(3.11) 



whereas £ is the restriction of £{t; z, 9) := K{t, r; z) dg^'^{T; z, 6) dr to (z, 9) = (1, 0). 
Using (|3.6p . also £ splits into pieces: £{t] z, 6) = X^Lo Z]p=o I^j=i ^^^iji't'i ^) ^fji^^ where 

/oo 
(-rf^-PKij{T;z)dg^''{T;z,9) dr. 
-oo 

We shall see that these t-independent factors are exponentially small with respect to g, as the 
latter tends to zero. Furthermore, denoting Kf-{t;z) := tPKij{t; z), 



1 i 2 



i=0 p=0 j=l 

This is so because, by Proposition 13.51 all integrals due to (13. lip factorize formally: if /i is a 
trigonometric polynomial at each order, then 



1 i 2 



poo '' ^ rco 

f {h£){T;z,9)dT = Y,Y.Y. f KT-z,9)TPK,,iT;z)dT 



i=Q p=0 j = l 



By virtue of Corollary 13.1 H we infer 

Proposition 3.12. On the homoclinic trajectory, i.e., setting {z,9) = (1,0), the following 
asymptotic expansion of the splitting matrix holds: 



/oo 
Xdef;{T)dT + Xc'^^{t)£f^, 
-oo 

where repeated indices are contracted (i S {0, 1}, p € {0, i}, j £ {1,2}J and 

poo p -. 

4j{t):=f [fl,^Y.^BM^^^K^{T)dT. 



(3.12) 



1=0 



As already pointed out, the terms appearing in the asymptotic expansion of Proposition [3?T2] 
will turn out exponentially small with respect to g. The factors shall pose no problems, 
the functions and K^^ being explicit and simple. In brief. Theorem H] begins to emerge! 

For the record, writing B^h{t;z,9) := j^_^K{t,T\ z)h{T; z,9) dr , we get (at {z,9) = (1,0)) 

/-t °o 
f:;,^Y.(^-^n'Kf^ {r)dT. (3.13) 

1=0 
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3.4. Emergence of exponential smallness. The "asymptotic" integrals in (j3.12p . including 
A Sfj , are of the form 

/OO °^ POO 

deF{X'';t-z,e)dt = ^e^ iqe^'^-^j [F{X'';t; z, q)Y dt, 

i=l 0<\q\<eN 

evaluated at {z,6) = (1,0), where the integrand on the right-hand side is a trigonometric 
polynomial of degree < IN by Theorem [21 We also point out that, by construction, the latter 
are ^-gradients, which allows us to omit the harmful g = terms. Above, F depends on 
X = X" locally, i.e., only through X{t;z,9), as well as analytically near X^, i.e., the series 



F{X) = F{X^ + X)=Y2 P^^^ {^Y 



(3.14) 



A:=0 



converges. In fact, since F is one of the functions in 

i^fl) {ji-K-tt'^^Kij^r I < p < i < 1 and 1 < J < 2}, (3.15) 

we observe that 

F(X«;t;z,0) =tPF(X«;0;ze^*,^ + c^t) (pG{0,l}), (3.16) 

with (z,6l) ^ [F(X";0;2,6')]^ analytic on (U^^^, \ {0}) x {|9m6l[ < a} for ah £ by virtue of 
Theorem [21 At z = there is a simple pole, due to -^02, at worst. 

We use the following lemma for analyzing such integrals. Its proof is given in Appendix lAl 

Lemma 3.13 (Shift of contour). Suppose that the function h{t;z,9) = tPh{ze^^,9 + tot) is 
analytic with respect to {z,6) G (U,-.^ \ {0}) x {|9m0| < a} and h{ - ,0) =0. Moreover, set 
tq := sgn(ci; • q)'dg~^ and 



o{n-<^~R)\t + ii,)P/i(e^(*+**'), g) dt 



(3.17) 



for each qe'L'^\{Q} . If 



zii-^\t + itqYh{e3^^+'^''\q)dt 



< A{g)e-'^\'^\ 



and if Hq{R) admits an analytic continuation to {0 < \R\ < p} with a pole of order k at 
R = 0, respecting the bound 

snp\Hq{R)\<B{g)e-''\'^\, 
\R\=P 

then we obtain the exponentially small (c > 0) bound 

k 



/OO 
h{t;l,0) dt 
-00 



< C 



With the aid of Lemma [3.131 we shall prove in Section [H the following key result: 

Proposition 3.14 (Convergence vs. exponential smallness). Fix at gM.. There exist positive 
constants c, ei and C , such that the estimates 

\T\t)\<cnJ- 

both hold true for all i = 1,2, .. .. 
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This dichotomy, in which the exponential smallness competes with the usual bound due to 
convergence, is not new; see [Gal94,GGM99b,Pro03]. Remarkably, we get the same exponent 
of the factorial, 4(z^ + 1)) as the latter articles — even though our method is quite different. 

Theorem [3] follows immediately from Proposition 13.14] bv an argument due to Gallavotti, et 
al: For each g, let n{g) be a positive integer. If |e| < ^ and e = eein{g)~^^'^~^^\ 

n{g) I I „ cx) I I „ 

\T{t)\ < Ce-'^a-'''-'^'' ^''^'^^'^ +C X; 

1=1 £=„(g) + l 

ei Vei/ 

since l\ < < n{gY in the first sum. Now, set n{g) = cg^'^^^^^^^ /In 2, such that n{g) ln|e| < 
— cg"^/'-'^^^-* and e = {e/ ei){c/ \n2)^'^'^^'^'> g^'^ (demanding in particular that eg~'^ be small). 

Remark 3.15. We used the exponentially small but diverging estimate to bound the partial 
sum e^T'"^) whereas the remainder of the series was easily controlled by convergence. As 

n{g) — > OO with 17 —s- 0, the important thing here is to have the exponentially small bound on 
for arbitrarily large £, in addition to the e-analyticity of T. 



4. Proof of Theorem [H 

We are left with proving Proposition 13.141 since Theorem H] was already shown to be its 
corollary. Here things are most conveniently explained using tree diagrams. However, each 
tree will be treated as an individual, solely for bookkeeping benefits, and no cancellations nor 
regroupings are forced upon them. 

By Proposition 13.12] we need to consider the simple factors c^j, as well as the integrals 

/OO °° POO 
t, z, q)dt^Y.''t [^(^"; ^' «)] dt (4.1) 
-co £ -j^ J — CO 

of Subsection 13.41 at each order i > \q\/N > 0. As F comes from the collection in (I3.15p . all 
integrals of the latter type shall be controlled with the aid of Lemma 13.131 

Due to the superscript u — referring to the unstable manifold — in the integrand above, the 
required bounds on the integrals over ]R_ in — f—00 ~^fo '^^^ straightforward, and are 
discussed later. 

In order to deal with we present the procedure below, which amounts to little more 
than integration by parts. First, we expand F according to ()3.14p and, like in Section [21 split 

X = X<i+ 62X = . + O ' 

dropping the superscript u from the notation. We can then express F pictorially as 

CO 

E E 

m=0 0<m'<m 
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where the binomial coefficient comes from the combinatorics of shuffling the arguments of 
the symmetric F^™) 0, which is attached to the root (node). Next, we replace the subtrees 

O = '^2^ with the expansion (j2.15p derived from (j2.1ip . 

In a tree T„g with root vq, a generic subtree T„ with root v has then the expression 

= U,{IC-'T^„. . . ; (X<i)^(""-'"")), (4.2) 

where T^. is a subtree — with root wj — entering v and the node function 

{jp{mv) cQj]2ing from (|3.15|) . if f = fo, 
y^(m„) ^ (j2.8p . otherwise. 

In a pictorial representation, there are lines entering the node v and precisely — m'^ of 
the latter are leaving a black dot (•) end node: 




The "whole" tree Ty^ contributes at orders i satisfying 

i>l + {m,, -m'J+Y, degT^,. =: d{T,,), (4.3) 

i=i 

where degT^, — defined in (12.140 — counts end nodes with suitable weights as well as nodes 
with exactly one entering line in the subtree T^, and the constant term 1 counts the root 
^jp{mvQ) — 0{e)). By this we mean that d{Tyg) is the largest integer such that 

T^^ = ©(edlT^o)) as e ^ 0. 

4.1. Some combinatorics. 

Lemma 4.1. The number of nodes, n{Tyg), of a tree Ty^ contributing at order i >2 obeys 

n(T,J<2(d(r,„)-l) < 2(^-1). (4.4) 
The number of end nodes is at most diTy^) — 1 < £ — 1. 

Proof. n(Tt,y) attains its maximum with respect to the "degree" d^T^^) as follows: 

(1) If n(Tyg) is even, there is only one line entering the root vq and the rest of the tree is 
binary, i.e., contains only end nodes and nodes with exactly two entering lines. 

(2) If n{Tyg) is odd, the tree is binary. 

Moreover, each of the end nodes is either • or (T), which contribute the least to the degree; see 
()2.14p . These choices minimize the number of end nodes when the 0{e) nodes having exactly 
one entering line are excluded (except at vq which is always 0{e)). Therefore, d{Tyg) gets 
minimized with respect to the number of all nodes, n(Tyf^). Since in a binary tree of j end 
nodes there are 2j — 1 nodes, we infer ()4.4p . 

The root and each of the end nodes increases d(r„g) by at least one, which implies the 
bound on the number of end nodes. □ 



This convention merely facilitates drawing. 
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Corollary 4.2. At most 2^^ trees contribute at order £. 

Proof. It is well-known that the number of (rooted) trees with k indistinguishable nodes is 
^(k) '■= ^(^k-i) — i^*^"^' '^tiich follows from generating functions (see [Drm04]) and the 
bound (2m)! < 4"^(m!)^. By Lemma l4.H the number of end nodes is less than £. We label 
the latter arbitrarily by the labels in {"lUK^ I 1 < < ^ — l| in order to form an upper 
bound on the number of our trees; these are the only possible labels, as otherwise d(T^p) 
certainly exceeds the order i — contradicting (j4.3p . The labeling of a tree with j end nodes can 
be carried out in at most (J'^^^) < 2^~^^~^ ways. The desired bound is thus obtained from 

1 + Efc=2 ^(^)(^fc-2^) + Cel2) EJ=7^^ ^(fc) < 2^^, because by LemmaOthe number of end 
nodes is at most i — 1 even though the total number of nodes can be as large as 2(i — 1). The 
term 1 on the left-hand side counts the single node tree F^^\ □ 

4.2. Simplification of integrals: scalar trees. We take a preliminary step towards bound- 
ing the values of the trees. 

Let us split the kernels K of the operators JC^^ appearing in (j4.2p into four pieces according 
to (j3.6p . These operators are attached to the lines between the nodes of a tree. Each of the 
2^^ trees counted in Corollarv 14.21 thus breaks into at most 4^^ new trees, as there are no more 
than 21 such lines by Lemma |4.1[ 

At the same time, we also expand the matrix products due to the coordinate representation 
of (1121) at each node: {U,y{fi, . . . , f„) = Ej+ii • • • Et=i(^-yn-jm^(' ' ' ' ^^^'^^ 
Tfc's represent all arguments of {i.e., subtrees entering the node v) — including X<i — and 
the superindices specify vector components. We next separate each scalar term into its own 
tree, thus getting up to {d + 1)^^ of these scalar trees from each old tree. 

There are lines carrying a factor t—r, coming from the i = 1 terms of p.6p . They correspond 
to double integrals: jJ_^{dT)'^ = dr (t — r). We remove these factors, insert a new node 
V on the line with a node function Uy = 1 and an integral sign both on the line leaving and 
on the line entering v. This operation can be depicted as 

^ 1 . (4.5) 

rrt f* r 

JJ—oo J— CO J— oo 

It does not affect the number of trees, but slightly simplifies the discussion below, even though 
the number of lines in a single tree can be as much as doubled. 
Altogether, we arrive at the following conclusion: 

Proposition 4.3. There are less than scalar trees to be considered, with = 2^^{d+ 1)^, 
at each order I. Each of these trees has at most 4{i — 1) lines. 

Remark 4.4 (Some conventions) . From now on, by a tree we will always refer to a scalar tree, 
where all the decompositions above have been carried out. In order to alleviate notation, we 
systematically omit all vector component indices. 

Recall that the node w' is the unique predecessor of the node w. Let Tyg be a generic (scalar) 
tree with root vq. Denote by V the set of all nodes and by T^nt the set of integrated nodes, i.e., 
the nodes whose leaving line carries an integral. We consider the root vq an integrated node, 
such that Vint consists of all nodes of Ty^ except black dot (•) end nodes. We can describe Ty^ 
by giving its structure recursively: if Ty is the subtree of T with root v E Vi^t, then 

Ty{t;z,9) = uy{t;z,e) J| / Ty,{T; z,9) dr. (4.6) 

w'=v 
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We call tit, the multiplier of the node v. It is a scalar function comprising all factors in the 
expression of the tree carrying the same variable of integration ("time label"), constricted 

in-between integral signs. In particular, it contains all subtrees • entering v, as these 

are the functions X<i involving no integrals. 

To be completely explicit, the multiplier of a node v € Vint is one of the functions below: 

(1) Kj2h^''''\ if V is the end node Q; see (l2lll . 

(2) Kyw'^'^^^ ( riw^Vint ) ( n^jeVint ) , if v is neither an end node nor the root vq^ and 



didn't appear from splitting a double integral according to (j4.5p . 

(3) 1, if V appeared from splitting a double integral; see (|4.5|) . 

(4) F(™") (n^^M„t ^<i) f n^-ev,, kJ) ,\iv = vo; see ^B- 



The functions and refer to diagonal elements of Kij and Kij in ()3.6p , respectively, m'^ 
is the cardinality of {w G Vint \w' = v} and m^j — m'^ the cardinality of {w G V\ Vint \ w' = v}. 
We draw the scalar (sub)tree in ()4.6p — originated from ()4.2p — as 




This diagram is reminiscent of the one below ()4.2|) . except that the operator has changed 

into the multiplier and the subtrees • have been absorbed into u^. Moreover, 

recalling the decompositions above, the present diagram carries on its lines instead of 
/C~^ — compare (j4.6p with (|4.2p — and has possibly more nodes due to the diagram in (j4.5p . 

We do not consider the end nodes with a black dot (•) nodes anymore. Subsequently, we 
will refer to "integrated nodes" (elements of Vint) as just "nodes". 

4.3. Integration by parts: one step. Reinserting the implicit arguments {z,6), 

u^{t; z, e) = tPu^{ze"'\ 9 + ut) (4.7) 

for the obviously defined Uv{z, 6). The power p can be nonzero only at the root, v = vq, where 
it possibly assumes the value 1 due to case (jH) in the list of all possible multipliers above. 
We now turn our attention to the Fourier transformed integrals 

rco roo rt 

1 Qvo ) I I / '^w 

{T;z,qw)dTdt, (4.8) 

l^'O+^^'w' =vo <1^=<1 w'=Vo 

which — recalling (|4.6p — arise from (|4.ip . at each order £ > \q\/N > 0. We omit the usual 
hat in Fourier transforms, since there is no danger of confusion. Due to Theorem [21 we may 
impose the finiteness condition 



\Qvo\+ 2^ <iN, (4.9) 

•w'=vo 

such that Ty^^{t; z, q) dt agrees at order I with 

'^es — / e^^'^ Uyo{t;z,qyo) TT / Tu,{t; z^q^) dr dt. (4.10) 

Here the asterisk reminds us that (j4.9p is being respected by the sum. 
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Remark 4.5. The task is to show that the t-mtegral in (j4.10p extends analytically from large 
positive values of 5ftc i? to a punctured neighbourhood of i? = 0, such that the residue can be 
computed. For this, we need the specific structure of the multipliers Uy. 

Let us define ^oi = Co2 = 1, Coi = Co2 = Cii = Cii = -1) and ^12 = ^12 = 0. There exist 
positive, continuous, functions aij and (iij such that 

\Kij{t;z)\ <a^j{z)e^'^''\^^ and \Kij{t; z)\ < aij{z)e^'^"'\^\ (4.11) 
hold on {(t, z) G M X C I |z| = 1 and ze'^* ^ ±i}. Moreover, 

aijaij < (4.12) 

where the constant A is independent of 7. Notice that, in the Kronecker delta notation, 

^ij+^ij = -2 5a5ji<0. (4.13) 

Emphasizing the node in question, instead of specifying the subindices ij we write and 
For each node v, we define the numbers 

ry:= sup mm{k eZ\y6 > : Uy{t;z,e)e~^^^^'^^^ ^ as t ^ 00} 
{z,9)eB 

with 

B := {\z\ = 1, |argz| < ^} x {|9me| < a}. (4.14) 

Thus, recalling (j4.7p . we are inside the analyticity domain U,-^^ x {|Qm0| < o"} of Lemma |2.2[ 
These numbers measure the divergence rate of the multipliers in the limit t — > 00. 

Lemma 4.6. Ordered according to the list of possible Uy 's on p. [221 

iv + kv, case (1), 

iv + {my - m;) + Y.w'=v case (2), 

0, case (3), 

Jf + {my - my) + Y.w'=v ^w, case (4), 

where S {0, 1} depends on the choice of F in ()3.15p . Moreover Uy{z, 6) , see ()4.7p . is analytic 
on {ISjm^l < 0"} with respect to 6 and on Vr,{f \ {0} with respect to z. It is also analytic in 
the punctured neighbourhood {\z\ > t^^} of z = 00, at which point there is a (possible) pole of 
order ry (if ry > 0). 

Proof. The maps ^^{z) and f{^^{z),9) are analytic in these domains, without singularities 
at z = 0, 00; see (jl.Sp . The rest follows by staring at the expression of Uy in each case. □ 

Starting from the end nodes of a tree — setting Sy ■= for them — we recursively define 

Sy ■= ^ ^ max(0, Ti^ ) and ny:=ry + Sy. (4-15) 

w'=v 

These numbers measure the divergence rate of (sub)trees in the limit t — > 00: 
Lemma 4.7. For Ty as in (|4.6p . and any 6 > 0, the estimates 

Ty{t; z, 0)e-(""+^)^* ^0 as t ^ 00, (4.16) 

(n r Ty,{T;z,e)dT^e-^'^+^^^^ ^0 as t ^ 00 (4.17) 

hold true in the region B of (I4.14p . 



ASYMPTOTIC EXPANSION OF THE HOMOCLINIC SPLITTING MATRIX 23 

Proof. Assume that (|4.16p is true for each successor w oiv {w' = v). By ()4.6p . we only need to 
observe that, given 6 > 0, \f_^ T^{t; z, 9) dr] < Ci + C2 /q e("»+^)T^ dr = C)(e["^^^(o,n™)+5]7i) 
in the hmit t — > 00. For then (j4.15p imphes (j4.17p for the node v. Now, (j4.16p follows from the 
definition of r^. Since (j4.16p holds for end nodes (n^ = 'Tv)-, induction proves the claim. □ 

Recalling Remark 14.51 l^t us now come back to the integral in (|4.10p . setting v = vq, i.e., 

e-^^Uy{t;z,q^) [J / T^{t; z, q^) dr dt. (4.18) 

An obvious problem is the exponential divergence of the integrand; see (|4.16p and (j4.17p . Our 
cure is the following. Since Uy is meromorphic at ^ = 00, by Lemma 14.61 we may expand 



Uy{z,9)= ^ z''uy^k{0) + uy,<^s{z,e) (4.19) 



k=~s 



for any integer s > —r^, with Uv,<-s{ztO) = 0{z ^ ) as z — cxd. Extending (14. 7p . we write 

Uv,k{t; 0) = e^"^* Uv,k{0 + i^t) and n^,^<_s(i; z, 0) = n„,<-s(2:e^*, 9 + ujt), 

where p G {0, 1} depends on the choice of F in case ([!]) on p. [211 In particular, the integral 
e~^* Uy^<,^sAt; z,qy)\[.^,^^ J_^Tw{t; z,qu,) dT dt is convergent for Kei? > -7, by virtue 
of (I4.17p . and can be estimated on a circle = p < 7 for the purpose of Lemma 13.131 

The rest of the integral (14.180 is integrated by parts: for —Sv<k<ry and sufficiently large 
positive values of 3?e R, 

foo rt 

/ e~^^ z^Uy^k{t;qv) W / Tyj{T;z,qyj)dT dt 

(4.20) 



2-^ oP-p' L 



\ p^^R^o-" Jo dt 



■dt 



where := R- k-j -iuj ■ q^ and Ey{t) : = Ey{t;z, {q^ \ w' = v}) : = n«,'=^ t-ao^^wir; z,qw) dr. 
But ^ = Y.w'=v'^w{t)Y{y,,=y^ij,-ty, f_oo'^w{T) dr , so that the recursion relation (jM]) yields 

^ = X] X] Uu,{t]z,qu,)(W I Tu,{T]z,qu,)dT\(^ j Tu,{t; z,qu,) drj . 

w-"" qw+Y.u,'=wi^=i^ «)'=«) w'=v 

We now fix w, and collect all the Fourier sums together so that they run over the set specified 
by qv + qw + Yjw'=w to + llw'=v,WTtw to = 1- Again, without affecting the £th order, we also 
impose the finiteness condition \qy \ + \qu,\ + '}2w>=w Itol + Yjw'=v,wi^w Itol < similar to 
(j4.9p . which allows us to bring the sum out of the t- integral. We observe that the remaining 
integral in (j4.20p produces integrals similar to the original (j4.18p . with the following changes: 
(i) Uy changes to u^, (ii) R changes to i?^, (iii) p changes to p' , and (iv) the integral on the 
line leaving w, J^_^Tu,{t) dT , is replaced by the ones on the lines entering w, the product 

Thus, a single step in the integration-by-parts scheme can be described in terms of trees as 
follows: Given a tree, consider one of the successors, w, of the root, v = w'. The line from w 
to V is "contracted" by erasing the node w, reattaching to v all subtrees originally entering w, 
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and replacing the multiplier of the root, v, by n^. Finally, we rename the root w. Pictorially, 





4.4. Integration by parts: exhausting the entire tree. We can proceed recursively and 
repeat the procedure in the previous subsection until there are no nodes left in the tree. Start 
with a tree having root v. First, call Tq := T„ and set fo := f . Then choose a successor w 
of V and define vi : = w, contracting the line from w to v. Next, in the new tree (the rightmost 
diagram above) called Ti, choose a successor of w and call it V2- Contract the line from V2 
to vi. Repeat until the tree has been exhausted and all nodes have been numbered. We can 
express the sequence of trees formed as 

Ti{t) = u,, Yl [ T^{T)dT, (4.21) 

wdTi: w'=Vi 



— oo 



where in the product we consider the tree Tj with the root Vi having the successors w. 
Let us define the numbers, the divergence rates of the multipliers tt^. , 

Ti := r„.. 

Analogously to the numbers in (I4.15|) . we set 

Si := max(0, n^). 

wdTi: w'=Vi 

The numbers are the ones defined in ()4.15p for the original tree, Tq. Notice that, although 
•So = , Si is not simply equal to , but is the analogue in the tree Tj of which Vi is the root. 
Similarly to ()4.17p . Si bounds the divergence rate of the product in (14.210 : if (5 > 0, 

H t r^(r)dr')e-(^»+'^)^* ^0 as t ^ oo. (4.22) 

The following, elementary, property is useful. 
Lemma 4.8. For every i, Sj+i + rj+i < Sj. 

Proof. Since r^+i := r^^^^, n^^^^ := r^^^^ + s^^^^, and T.weTi:w'=v,+i max(0,n^) = Sy^^^, and 
Sj-|_i := max(0,n^)= max(0,n^)+ max(0,n^), 

we get Sj+i + rj+i = Y.weTi:w^Vi+i,w'=Vi (0' + '^"i+i - ^i- Here we used the fact that 
v'-_^_^ = Vi in the tree Tj. □ 

Let A;o : = and Qq := 0. Suppose that at the ith. step we are considering the integral 

f'OO _ rt 

/ e(*«-'^+'='->-^)*t^'»n,,(ze^*,g,Je*«"---* J] / T^{t; z,q^) dr dt; 

Jo ^rr I J —OC 

of. (j4.18p . Then the integration-by-parts procedure described above takes us at the (i + l)st 
step to a similar integral with i replaced by z + 1, defining 

ki+i ■.= ki + ki and Qj+i ■ = Qi + Qvi, (4.23) 
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multiplied by the factor 



[R - ki+ij - iLO ■ Qi+i) 



(4.24) 



The integer indices Pi+i and ki can assume the values 

< Pi+i < Pi < 1 and ki = -ki - Si,..., Vi. (4.25) 

With the aid of Lemma 14.81 it is straightforward to see that 

< A;i + r, + Si < n^o, (4.26) 

where n^^ = s„g + r„Q = sq + tq is the number describing the divergence rate of the original 
tree, Tq, in the sense of (I4.16|] . The number of possible values of ki above is thus at most 
1 + TiyQ. This implies the two-sided (equivalent) bounds 

ri - n^o <h<ri and < - fcj < n^i^ . 

We continue recursively until there are no nodes — alternatively, integrals that require reg- 
ularization — left. The result is a sum with terms of two different species. We set 

Ri := R - kij - ioj ■ Qi, 

in order to make the presentation more compact, and also define (c/. E^j in (j4.20p ) 

Ej{t) := Ej{t;z,{quj\w £Tj,w' = Vj}) := / T^(t; z, g^,) dr, 

with the understanding that E^y.^^^_i = 1. The first class of terms is 

f n ^^^fe^l r -.,,<-h-sM^'^^^^E^^ dt, (4.27) 

for < j < |Vint|- Notice that s\y.^^^\_x = 0. Second, there are the terms (0 < j < |Mnt|) 

I ii pi+p»-pi+i I pi+Pj ' ^ ' 

\ i=0 / ^j+i 

4.5. Estimates. By the definition of Sj, we can bound 

\Ej{t)\ < Ci=;^.e(^^+3)^*e-'^^l9»l {t > 0), 
for \z\ = 1 with |argz| < i?. Deferring the proof to Appendix lAl we formulate 
Lemma 4.9. The coefficients Ut,.^^. (g) satisfy 

k.,,,,(g.J|<a,r'=-'^^e--l''"»l, 
where r is as in Lemma \4-6[ On {1 < \z\ < t^^ with [argz[ < i?}U{kl — 
K,^^-,^_,S^,q.J\ < a.(r|z|)-^'-^-i2"«or->e-'^l'^"«l. 
When \z\ = 1, |argz[ < ■!?, and |9m0| < a it holds true, for suitable fi G Z, that 



\u,,^{ze^'\6 + u;t)\ < 



e^^-y^ t < 0, 

e*^'^*, t > 0. 



26 MIKKO STENLUND 



Recalling K2E\\ . the integral in K27\\ is bounded by (2r-i)""o+i%^e"'''^"j ' 



7 '"^ 



when 3?ei2 > — A7, 1^1 = 1, and [argz| < ■!?. Thus, the bounds on (j4.27|) and (j4.28|) read 



2 



i=0 



and 



j=0 [-"'i+l I / 



respectively. Here e J^i^ol'Ji'J ^^El*?™! < g o"!?!^ because the tree is a Fourier transform with 
index q, and the g^. , come from convolutions. The factors of r are controlled by 

Lemma 4.10. For all j, X]i=o(^* ~ ^i) ^ —^vo- Moreover, Uy^ < d{T.u^^) (see ()4.3p ). 
Proof. From (j4.23p and (|4.25p . we clearly have fcj+i > — Sj. Then, 

y^(fe» - -Ti) = A^j+i - ^ > -Sj - - ^ rj > -Sj-i - ^ rj > • • • > -sq - ro = -n^g, 

i=0 i=0 j=0 i=0 

with the aid of Lemma l4.81 The estimate n^^ < d(r^J follows easily from g3]) and (ITO]) . □ 

By the last bound of Lemma 14.91 it is clear that Cej = (nl=j+i^ ^fO^^j' fact, each 

integral in Ej produces a 7"^ in the upper bound, and we get Ce^ < (C7^^)l^'"'l^''^^. Starting 

from the explicit expressions of the multipliers, we estimate l\i=o^ ^ < {Cg)^^^{Cef^^-o') 
with the aid of ()4.1ip and ()4.12p . Here V, defined above ()4.6p . is the set of all nodes including 
the black dot (•) end nodes. Since \Viat\ < \V\ < 2n(Ti,„) < 4(d(Ti,g) — l). 

Notice that the fraction 1/Ri is analytic on the domain where \R\ 7^ \ki'y + ico ■ Qi\. Thus, 
all the 1/Ri are analytic in the punctured neighbourhood of the origin 

< < p:= imin{{7}U{|L^-Qi| | Qi^O}}. 
On the circle \R\ = p, < p^^^^, provided p > —1. According to (I4.25p . 



n 



1 



pi+pi-pi+i 



<p'J-P0+Pj (0<i<|Mnt| 

1=0 

Subsequently, summing over j the bounds on ()4.29p and (I4.30p we get 

i=i 

At the ^th order, \Qi\ < £N for each i, and thus p > ^ min(7, a(£A^)~'^). 

From the above it is also clear that the order of the pole at i? = in our integrals does not 
exceed |Mnt| + 1- For the purposes of Lemma [3. 131 we compute (using 5,7 < C) 

j=o ■ ^ 
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There are at most (4^)! < 4^^(^!)'^ orders in which we can exhaust ah of the up to Ai hnes of 
a tree contributing at order i. 

4.6. Remaining integrals. The integral over M_ in (14. ip is simple, because the integrand 
satisfies the identity ()3.16p and has the analyticity properties stated below that equation. In 
other words, we can separate the (possible) pole and the constant term from the rest: writing 
FiX^;O,z,0)=F{z,e), 

Fix""; t, z, 6) = tPz-ie-T*F_i(0 + ujt) + t^FoiB + cot) + tP5iF{ze'\e + ujt). 
Here p = 0, 1. Applying on each term separately, we get that j-^^ F{X^; t, z, 9) equals 

These terms are small compared to to the large bounds obtained for the part above. 

We also have to study the coefficients c^j appearing in ()3.12p . This is most conveniently 
done in terms of the representation (|3.13p . Since = 0(e), only 



(T)dT 



with / < i can contribute to c^j at order i. We only need to consider t > 0, as T(t) = T(— t). 
The above integral consists in obvious shorthand notation, through (|3.6p . of 4' terms like 



£ 



{~t - nr^^ {fKi){riU in - T,_i)f' (i^,M)(rz_i) • • • Xi(ri)-f (n - Tof^ (i^i Af)(ro)TP"i^o(ro), 



Pa , 



where pi^i = 0. If = 1 (instead of 0), we use tq = (tq — ti) + ■ ■ ■ + (t;_i — r;) — (— t — r^) — t, 
getting / + 2 terms of the original form except that pq = and either there is a factor t or 
precisely one change pi i-^ pi + 1 occurs. Due to M = 0{eg^), (j4.1ip and (|4.12p . each KiKiM 
produces a factor CAIelf?^'"*"^ to the upper bound, whereas analyticity yields e"'^'"''! (which we 
prefer although at each order we reduce to trigonometric polynomials); ni=o — c"'^''^'. 

Last, the integral involves a total of ^'=}(1 + Pi) < 2/ + 3 denominators of the form 
ki'y + iuj ■ Qi, where ki € Z and Qi := J2]=o with Qi = q ^ 0. We bound these by 

ymm{{j}U {\uj ■ Qi\\ Qi 0}} J 

Altogether, the bounds above yield C^^\e\'-g-^e-''^'i^{i\f''. 

Skipping further details, this is the upper bound on the integral. It is smaller than what 
was derived for the integral £^ F — and thus for £fj in (j3.12p — above. In conclusion, among 

the contributions (4) (ff^) (h + £2 = i - 1) to {Sfj) is the most dan gerous one. 

Remark 4.11. Above, the sums over q^, with Y2i1i ~ Q — were dealt with as 

follows. Since the analyticity domain with respect to 9 is the compact {IQ'm^l < a}, it can be 
substituted by some {|9m0| < a'} with a < a'. Then, for every i we actually have the factor 
g-f^ M in the estimates above. While e~'^^»l'''l < e^'^l'^', we get rid of the sums over qi: 

A similar remark applies to the Fourier indices in Subsection 14.51 □ 
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5. Discussion 

Each term in the asymptotic expansion (j3.12p of the sphtting matrix T is proportional to an 
integral of the form £^ dgF^X"^; r, z, 9) dr. These we showed to be exponentially small in the 
limit 5 — > at all orders i, by extending the integrands analytically into a wedge {|arg2;| < i?} 
on the complex plane and then shifting the contour of integration. A key point is that we had 
to do this order by order, because the series 

E£oe^^"'^(^'^) is not expected to converge for 

large values of |z| if e is fixed. 

The large powers of the factorial i\, associated with the regularized integrals, are produced 
by accumulation of poles at the origin in the R plane. To some extent the factorials are 
artifact, as is shown by the following simple example. In order to study, say, the integral 

oo /-t 

u{e + uJt) / v{e + uJT)ze^^ dr dt, 

J — oo 

we have to show that the sums 

/•oo pt 

Vn(p)?}(g-p) / e^ip-^-^)t e^'^'i-P>'^+^^^ dT dt 

p Jo J-oo 

extend analytically to a (punctured) neighbourhood of the origin, R = 0. We integrate by 
parts, just as in ()4.20p when extending analytically the tree integral ()4.10p . and get 

y2u{p)v{q-p)-^ — - — \ — — I ■ — - + - — 

^ R-ip-uj \^{i{q-p)-uj+-f) R-ij + iq-u) 

Here the pole at ip ■ u! gets arbitrarily close to the origin, unless we restrict p somehow — for 
instance, by considering trigonometric polynomials. Either by simplification, or by computing 
the same expression directly by starting from the inner integral, we obtain 

1 — • - — -. 

^ i{q-p)-uj + -f R-{j + iq-u) 

In the latter form there is no problem; the pole has cancelled. Of course, this is a naive example 
and in general it is hard to see whether a given pole popping out of the integration-by-parts 
procedure should really be there. 

Also the coefficients cf^ appearing in Proposition 13.12) produced large powers of ^!. Even 
though integration by parts was not exploited, the source of the factorials was again the 
accumulation of poles at the origin in the R plane. In both this case and the previous, the 
type of "divergence" is very similar to what is encountered in KAM theory. There repeated 
resonances, or arbitrarily many occurrences of the operator T>~^ in convolutions, ruin absolute 
convergence of the Fourier-Taylor expansion of a solution by producing high powers of the 
factorial i\. On the other hand, the state of affairs can be cured by well-known resummations, 
as in [Gal94]. Such resummations still escape us in the context of homoclinic splitting. 

Appendix A. Some computations and proofs 

Proof of Theorem O It is known that the whiskers are Lagrangian manifolds [Eli94, DGOO, 
LMS03]; they are graphs over the angles and, moreover, the actions are gradients of potentials. 
In other words, there exist generating funcions S"'* such that 

r'^ = a<^5"'"(0,V') and A"'^ = a^5"''((/),V'). 
We consider the latter row vectors. Differentiating both sides of 
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with respect to the ith component of {(pjtp) =: (p = (930, ^i, ■ ■ ■ , ^d) yields 

d d 

j=0 j=0 

Here J := {I, A) and Qj^ W"'^ means Kt^^j = aj^7^(</., V, 9^5"'^). 
From (lA.ip it follows that, on the homoclinic trajectory, 

dj7e''dl{S'' -S') = Q (A.2) 

because one has 9(^5" = d^pS^ and therefore also Ti^ = for i = 0, . . . , 2(i + 1. For our 
particular Hamiltonian, djTC^'^ = (I, A)"'*, where we now drop the superscripts u,s and 
evaluate everything at the homoclinic point 

^ = (<^, ip) = (0, iot) + X"(e^*, wt) = (0, ujt) + X'{e"'\ujt) (A.3) 

below. The ■0 component of ()A.2p reads 

19^945" - 5^) + (5" - 5^) = 0, 
or, recalling / = /° + 0{e) = 2g/ coshgt + 0{e) ^ (taking e := g~^e^l*le small), 

5^5,^(5" - S') = -r^AdliS'' - S'). (A.4) 

Here A is the row vector {Ai, . . . , A^^), such that both sides of the equality are row vectors. 

Next, let us perform the coordinate transformations {(pjip) = F^'^{z,9) 1— > {z,9). To this 
end, we observe that the splitting vector satisfies 

Then, the 9 derivative of the column vector A is the square matrix 

by the chain rule. In particular, by (jA.4p . 

d^A = - S')f = -5^(5" - 5^) r'A^ = -d^A r'A^ (A.5) 

holds at a homoclinic point (see ()A.3j) ). such that 

deA = d^A{de^-r^A'^decP). 
The matrix M : = deij) — I'^A^decj) has the asymptotic expression 

M = (1 + 0{e)) - + 0{e)) • (a; + 0{e)f ■ 0{e) 

with e = 5e~^l*le, as before, and e small. □ 

Proof of Lemma \3.1l (Adapted from [SauOl]). The Fourier transform of ()3.2|) yields 

F{s, q) = F{0, q) e-^(9"'-'?)^, (A.6) 

which is entire in s. But |F(0, g)| e^f"''^ '?)^'"'' = \F{s,q)\ < B{g)e-'^\i\ for s G and 

\F{0,q)\ < B{g)e~'^3-'\^-g\-v\g\_ 

Finally, plugging this into (|A.6p . we get 

\F{s,q)\ < 5(5)e-(''-l^'"^l)3~'l-9l-^kl (s ^ C). 

For |9ms| < '0 and |Qm6'| < rj' < r], the series J2q£Z'' -^(^' q)^^'^^ is uniformly convergent and, 
as such, provides the analytic extension. 
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Since ax ^ + /3x > a{v + 1)(^) ^^'^^'^^^ for positive a, j^, and a:, we get from the 
Diophantine condition (|1.4|) that 

holds if s G M, g G Z'^ \ {0}, < 5 < r/, and i(;(??,7/ - S) := (??a)i/(^+i) (77 - 5)^/(^+1) (1/ + 
Moreover, ^gg^d^jo} e^'^l'^l < C6~'^, where C only depends on the dimension d. 
By (iX6l) . -F := • )) = F(s, 0) = F(0,0), such that F does not depend on {s,e). □ 

Proof of Lemma \3.1^A By shifting the contour of integration from R to the complex plane by 
itq := isgn(u; • q)i!}g~^ units, we compute 

1 



f 

J — ( 



tPh(e3\q)e''^-'^^dt= res „ 
R=0 R 



-•dg ^\Li>-q\ 



where Hq is defined in (ISTH) and Iq{R) := e(*<?-'^+^)*(t + itqYh{eS^^+'^'i\q) dt. There are, 
a priori, two additional line integrals JJ*'', but they cancel due to the residue at i? = 0, as is 
easily checked. Because /i( • , 0) = 0, Iq{R) does not have a pole at i? = 0. Hence, 

e'^^iIq{R) 



res 



R 



z'^-'^\t + itqYh{e3^*+'^^\q)dt. 



If Hq(R) has a pole of order k at R = 0, then 



res ■ 
R=o 



Hq{R) 



R 



E 

j=0 



i-itqyHq^-J 



<E4(^y-piF,(i?)i 



because the Laurent coefficients Hq^^j := ^ !f\R\=p ^^-^ can be bounded from above by 

\Hq-j\ < sup|^|=p|-ffg(i?)|, whenever the circle \R\ = p is inside the domain of Hq. 
Under the assumptions of the lemma, for any < 5 < a and q ^'L'^\ {0}, 

k 



r rh{e3\ q)e"^-^' dt < U^) + B{g) ^ ^ (^)' 



by mimicking the proof of Lemma l3.1i Here w{'d,a — S) := {'daY^^'^^^'' [^-^Y^^'^^^\iy + 1). 
Summation over q produces a factor C6~'^. □ 



Proof of Lemma \4-9[ Here we usually omit the subindices i,v, and vq. According to Lemma HTG 
is analytic on {\z\ > t~^} x {|Qm0| < a} with a pole of order r at z = 00. A Cauchy 

estimate then reads \u.v,k{d)\ < t'' sup^^^^^-i\uv{z,6)\ < Ct^~^. Now the first bound follows 

from \uy^k{q)\ < e"'^l^l sup|Q^e|<^|u^,fc(6')|. 

The second inequality is a trivial consequence of the first one for \z\ > 2r^^. For the other 

values of z, one uses ()4.19p and ()4.26p as well as |2:|~''|u^(2:, g)| < Ce"""'*^' to bound 



\u. 



v,< 



.k-siz,q)\ < \z\ 



-k—s 



i 



r+k+s 



< Ce-''\''\\z\ 



-fc-s-l 



(2r-^) 



i=0 

— l^fc+s+r+l _j_ ^— fc— s— r— l2'"+fe+s 



The last inequality in the lemma is obvious and is stated for completeness. 



□ 
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